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Abstract 

A fiber composite flat cantilever plate that 
has symmetric and nonsyrametric laminate configura- 
tions is theoretically investigated to determine 
its static and dynamic structural response. The 
finite element analysis method used includes a 
unique triangular finite element developed at Lewis 
for the analysis of fiber composite airfoils. The 
various responses investigated include tip dis- 
placements, natural frequencies, and fundamental 
mode shapes. The results show that laminate con- 
figurations may be selected for a cantilever such 
that when the tip at the leading edge is loaded 
normal to the plane of the plate, the tip at the 
trailing edge can (a) deflect in the opposite di- 
rection, (b) deflect about the same, or (c) deflect 
more than the tip at the leading edge. This varia- 
tion in response can be utilized tp provide built- 
in structural damping to resist flutter. The re- 
sults also show that the displacements and the 
natural frequencies can be in considerable error 
for nonsymmetric laminate configurations if the 
membrane -bending coupling is not taken into ac- 
count. Structural response results obtained from 
ten laminate configtirations are presented in tabu- 
lar forms and may be used as an aid in selecting 
laminate configurations for composite airfoils. 

Introduction 

Fiber/reain composite laminates are being used 
or proposed for use in structural components such 
as wings, horizontal and vertical stabilizers for 
airplanes, helicopter and wind-energy machine 
blades, and fan blades for aircraft turbine en- 
gines. Laminates for these types of components 
have a large ntimber of plies oriented at different 
directions and are selected to meet design require- 
ments for strength and stiffness. 

Current practice is to design the aforemen- 
tioned components using laminate configurations 
T^iich are balanced (no in-plane shear -stretch- 
coupling) and are symmetric (with respect to bend- 
ing) . One reason symmetric laminate configurations 
have been studied extensively for these components 
is that they are amenable to solution using well- 
known analysis methods ► »2) Oq the other hand, 
unbalanced and/or nonsymmetric laminate configura- 
tions may possess shear -stretch, bend-stretch , and/ 
or bend-twist coupling responses (Appen. 1 of 
Ref. 3) . Laminates with these coupling responses 
are difficult to analyze and partly for this 
reason, have been avoided for use in structural 
components. 

Some important reasons for investigating lami- 
nates having coupling responses are: 

1. It may be possible to select coupling in 
laminates that will provide built-in self -damping 
fflfic ban isms when subjected to dynamic excitations, 

2. Structural components fabricated from fiber 


composite laminates are occasionally found to be 
warped and it is believed that the warpage is 
largely due to coupling responses that result from 
fabrication errors in orienting the various 
plies. (5) 

To date, loading displacements of cantilevered 
laminates with coupling responses have not been in- 
vestigated theoretically. In addition, the vibra- 
tion responses of these laminates have received only 
limited theoretical investigation using primarily 
indirect methods.^^) 

The main objective of this investigation was to 
Investigate the displacements and vibrations of 
cantilevered flat laminates with various types of 
coupling responses. TWo secondary objectives were: 

1, To identify some coupling responses which 
may give rise to built-in self-damping mechanisms, 

2. To indicate how molds may be contoured to 
permit, or restrain, warpage resulting from fabri- 
cating laminates that have various co\jpling re- 
sponses . 

An analysis method was used \^ich utilizes a 
finite element developed at Lewis for the analysis 
of aircraft gas turbine fiber composite fan blades. 

A summary of the equations defining the element 
stiffness matrix is given in the appendix. 

The cantilever used in the investigation is a 
rectangular plate with aspect ratio (a/b) equal to 
2 (Fig. 1).. For the displacement response, the 
cantilever was loaded with a concentrated load nor- 
mal to the plane of the plate at the leading edge 
tip. For the vibration response, the cantilever was 
free of loads. Fiber composite laminates with vari- 
ous coupling responses were selected to illustrate 
the effects of these coupling responses on the dis- 
placements and vibration frequencies of a canti- 
levered plate. 

The effects of the coupling responses were also 
investigated using the reduced stiffness (reduced 
bending rigidities) concept (Ref, 3, Appendix 1). 

Laminate Geometry and Configuration 

The structural component chosen for this in- 
vestigation was a cantilever plate represented 
schematically In Figure 1. The reason for select- 
ing this component is that it conveniently simulates 
an airfoil. The aspect ratio (length /width) equals 
two and the width -to- thickness ratio equals ten. 

The dimensions of the cantilever axe a *= 2 inches, 
b = 1 inch, and h = O.L inch (Fig. 1), 

Ten different laminate configurations were se- 
lected for Investigation. The ply orientations for 
these laminate configurations are shown schemati- 
cally in Figure 2; twelve plies of equal thickness 
were chosen for these laminates. Fotir cases are 
ayiiiinetric with respect to bending (cases I to III, 
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and V); the other six cases are nonsyinmetric. Note 
that case I has all 12 plies at 0^ and case V has 
6 plies at 0^^, and 6 plies at 90°, These are the 
only cases with orthotropic material behavior. 

The composite material chosen was Thornel- 
75“S/epoxy with a fiber volume ratio of 0.6. The 
plate constitutive properties (force-deformation 
relationships At j , Ci j , and Dij) were computed 
using the program described in Reference 6 and are 
tabulated in Table 1. (Also see the Appendix.) 

The Aij , Cij , and Dij matrices are symmetric , 
thus, only the upper triangular portion is given. 
Note that for the symmetric laminates (symmetric 
with respect to bending) the Cij* a are zero, char- 
acteristic of the fact that no bend-stretch cou- 
pling exists for these cases. Note also that all 
cases except cases 1 and V in Table 1 have either 
some nonzero Cij terms, Ai 3 » A 23 terms, or Di 3 , 
D 23 nonzero terms. These cases exhibit anisotropic 
material behavior. In the reduced stiffness con- 
cept the Aij and Cij arrays are incorporated 
into Dij via the following: 

Dij ^ ^ij ” where ^ij the laminate 

reduced bending stiffness . The reduced bending 
stiffness coefficients for dJj are given In 
Table 2 for the laminates Investigated. 

Method of Analysis 

The analysis uses a unique triangular thln- 
plate finite element. This element is deduced 
from a thin-shell, double-curvature, variable- 
thickness, isoparametric, anisotropic finite ele- 
ment which was generated for the analysis of com- 
posite airfoils (Fig. 3) . The element has six 
nodes ( 3 -corner, 3 -midstde), each with five dis- 
placement degrees of freedom (DOF) , consisting of 
three translation displacements, u, v, and w, and 
two rotations, a and The displacements within 

the element are represented by parabolic interpola- 
tion functions. The coordinate system within the 
element uses the same Interpolation formula as the 
displacement; therefore the element is isopara- 
metric (see Appendix). The a train -displacement 
equations used in the formulation are from thin- 
shell theory and include double -curvature, bend- 
stretch coupling, and through -the -thickness shear. 
The assumption that plane sections remain plane 
after bending is relaxed. The equations for the 
element stiffness matrix are summarized in the 
Appendix . 

The thickness of the element ia an input 
parameter at each node. The material coupling re- 
sponses are represented by the twenty four coeffi- 
cients which are required to describe the force- 
deformation relationships within the element. 

These coefficients are elements in the Aij, Cfj, 
and Dij matrices mentioned previously and tabu- 
lated in Table 1 for the cases Investigated. The 
element stiffness and master (global) stiffness 
matrices are assembled In the conventional way to 
yield the global system of equations F - Ku, 
where F represents the nodal forces, K the 
master stiffness matrix, and u the nodal dis- 
placements. The solution of the system F = Ku is 
obtained by solving this system for the displace- 
ment variables (u) . 

The consistent mass matrix (M) is generated 
using the previously mentioned parabolic interpo- 
lation functions. In addition the mass density of 


the element can be a variable and Is therefore given 
as input at the individual nodes. The matrix equa- 
.2 

tion, Ku + M — j (u) >" 0, is solved for the elgen- 
dt^ 

values and eigenvectors using available eigenvalue 
extraction routines. This allows the calculation 
of the natural frequencies and the normalized mode 
shapes for structures that the finite element has 
been intended to represent. The element formula- 
tion was checked against known solutions for both 
displacements and natural frequencies. 

Results and Discussion 

The results discussed herein include nodal 
diaplacements due to a concentrated load normal to 
the plane of the plate at the leading edge tip, 
using both 5 and 3 DOF models with 4 and 16 ele- 
ments (Fig. 4) and also the first six fundamental 
frequencies (harmonics) with their corresponding 
mode shapes (Fig, 5). Results from laminates with 
coupling are compared to corresponding results ob- 
tained using the reduced stiffness concept. 

Catllever Tip Deflections 

. Table 3(a) gives In tabular form the tip de- 
flection from the leading edge (LE) to the trailing 
edge (TE) at the quarter stations (Fig. 4(b)) for 
the cantilever (Fig. 1) without bend -stretch cou- 
pling (only bending variables) , This model con- 
sisted of 16 elements (Fig. 4(b)) with three dis- 
placement DOF (w, a, 3) per node resulting in 120 
free variables. The applied load for all the dis- 
placement results was a 15-pound concentrated load 
at the tip at the leading edge . The aluminum case 
is Included to illustrate the corresponding be- 
havior of a typical isotropic material. Table 3(b) 
gives the displacement results with five displace- 
ment DOF per node (200 variables total) • The affect 
of bend-stretch coupling (nonsymmetry) is evident 
when the displacements for cases IV, VI to X given 
in Table 3(b) are compared with the corresponding 
results given in Table 3(a) . It was previously 
pointed out that case I is the orthotropic case 
where all 12 plies are in the longitudinal direc- 
tion. This Is the stlffest case and exhibits a 
slight downward displacement at the trailing edge, 
Case X has very slight material bend-stretch cou- 
pling but it is noted that the trailing edge has a 
displacement larger than that obtained at the lead- 
ing edge where the load Is applied. With the ply 
orientation of cases IV and X, deflections are ob- 
tained which are very large compared to the ortho - 
tropic case. Table 3(c) gives the displacement re- 
sults for the material coupling cases using the re- 
duced stiffness values in Table 2. Only the bend- 
ing variables are used with the reduced stiffness. 
Thus, the total number of variables decreases from 
200 to 120. As can be seen the results obtained 
are comparable to those given in Table 3(b). 

The previous discussion leads to the following 
important observation. Since the deformed shape 
for unsymmetric laminates is calculable, laminate 
configurations may be selected with predetermined 
twist (untwist) for anticipated membrane and bend- 
ing loads . This can be used to offset Increasing 
angle of attack in airfoil designs and thereby pro- 
vide structural damping to minimize or avoid 
flutter. Also the predetermined deformation can be 
used to contour the laminate fabrication mold. 
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Cantilever Frequencies 


Table 4 gives the first six natural frequen- 
cies of the cantilever plate for the ten different 
ply orientations. Table 4{a) presents the results 
using only bending variables, w, ct, 0 with the 16 
element model and 3 DOF per node. For Table 4(b), 
the element model is reduced to 4 elements 
(Fig, 4(a)). In Table 4(c), results are given for 
a 4 element model with 5 DOF per node. Comparing 
corresponding results in Tables 4(a) and (b) , it is 
seen that there is very little difference between 
the results using the 16 element and the 4 element 
models for the first and second natural frequen- 
cies . The 16 element model usually predicts lower 
frequencies than the 4 element model for the cases 
without coupling. For the higher frequencies, 
there is a larger difference with the 4 element 
model generally underestimating the natural fre- 
quency* Note that the natural frequency of the 
ortho tropic case is the highest for first frequency 
(first bending) but is low compared to the other 
cases for the second frequency (first torsion) . 

The fifth and sixth frequencies for cases II to IX 
are higher than those for case 1, For the cases 
with material coupling (Table 4(c)) all the modes . 
are lower than those without coupling (Table 4(b)). 
The maximum difference that ia obtained for the 
first frequency (using the case with coupling as a 
baseline) is 72 percent in case VI, while there is 
a 91 percent difference in the sixth frequency for 
case X. 

Cantilever Frequencies Comparisons 

A comparison of the vibration frequencies of 
the cantilever plate with and without coupling and 
with reduced bending stiffness is given in Table 5 
for the 16 and 4 element models for case VII, For 
the 5 DOF per node idealization, the 4 element 
model first frequency is only 1,5 percent higher 
than the 16 element model. The largest difference 
is 13 percent and is in the third frequency. For 
the 3 DOF per node idealization (no material cou- 
pling) both the 16 element and the 4 element models 
estimate the first frequency to be 44 and 46 per- 
cent higher, respectively, compared to- the 5 DOF 
Idealization. All the results obtained using the 
3 DOF idealization without material coupling over- 
estimate the natural frequencies with a maximum 
difference of 60 percent for the sixth frequency 
compared to the 5 DOF idealization. The 3 DOF per 
node idealization, therefore, is not a realistic 
representation for a laminate having material cou- 
pling, This idealization does not have the 
stretching (membrane) flexibility that the u and 
v variables provide the formulation. Therefore, 
membrane flexibility is necessary to Insure accu- 
rate frequency results . 

The results presented in column 5, Table 5 
are for the 3 DOF per node, 4 element model using 
the reduced stiffness values for the material con- 
stants in Table^2. The first and second frequen- 
cies are only 3 percent and 7 percent lower, re- 
spectively, than those of the 5 DOF per node 
16 element model. The sixth frequency differs by 
less than 6 percent for these two cases. The re- 
sults indicate more flexibility for the first two 
frequencies and greater stiffness for the higher 
ones. Case VII was chosen for these comparisons 
because it had significant material coupling as 
shown in Table 1. 


Table 6 compares frequency results with and 
without simulated centrifugal stiffening for the 
5 DOF per node Idealization, 4 element model. Sim- 
ulated centrifugal stiffening was introduced by 
setting the displacement variable u In the x- 
direction (Fig, 4) equal to zero at the cantilever 
tip. Note that the first frequency Is increased by 
20 percent. The remainder are increased by 7 per- 
cent or less with the exception of the fifth fre- 
quency which Is Increased by 24 percent. Compari- 
son of these results show that centrifugal stiffen- 
ing affects the first frequency but has little 
effect on the next three higher frequencies. 

Cantilever Vibration Mode Shapes 

A pictorial representation of the first six 
mode shapes is shown in Figure 3. The mode shapes 
are normalized with respect to the leading edge tip 
displacement. Note that for case I, the third mode 
is a transverse bending mode which does not appear 
in the isotropic material (aluminum) until mode 6. 
Also note that the third bending mode of the iso- 
tropic material does not appear in the first six 
modes of case I. Case VII was chosen because it 
has sigriificant material coupling while case IX has 
strong bending stiffness and complete coupling 
characteristics. For case VII, modes 4 and 5 are 
.very close numerically (within 6 percent) as shown 
in Table 5 for the 5 DOF per node, 16 element 
idealization. The pictorial representation for 
these nodes is quite similar . Compared to the iso- 
tropic material, there is a slight shift in the 
nodal line for the second bending mode; the third 
bending mode appears as the sixth mode for cases 
VII and IX. 

Cantilever Tip Deflection Comparisons 

Table 7 compares the tip displacements ob- 
tained for case VII using both membrane coupling 
and reduced stiffness with the displacement ob- 
tained without coupling (3 DOF per node idealiza- 
tion) . As can be seen, the differences between the 
coupled and uncoupled cases are from 42 percent at 
the leading edge to 55 percent at the trailing edge. 
Using the reduced stiffness values for the material 
constants, the variation is only 6 to 7 percent. 

This shows the artificial stiffening that is de- 
veloped if the idealization is not permitted mem- 
brane flexibility. The reduced stiffness is shown 
to give very good displacement results compared to 
the 5 DOF model, which is believed to be an ade- 
quate representation of the cantilever studied. 

The previous discussion leads to the following 
observation. Use of the reduced stiffness concept 
to account for bend -stretch coupling in cantilevers 
made from unsymmetrlc laminates provides a good 
approximation for predicting displacement and vi- 
bration frequencies. Therefore, finite element 
computer codes without bend -stretch coupling can be 
used to analyze nonsymmetric laminates using the 
reduced stiffness approach. Compared to the case 
where the bend -stretch coupling is included in the 
formulation, the reduced stiffness approach has the 
following advantages: (1) Uses simpler finite. ele- 

ment formulation, (2) requires only six stiffness 
coefficients In the force-deformation relation- 
ships, (3) For the same finite element representa- 
tion of a structure, it will take about one -third 
the computer storage and will run approximately 
five times faster. 
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Summary of Results aad Conclusions 


The major results and conclusions obtained 
from this investigation dealing with the structural 
response of unsyinmetric anisotropic, flat, com- 
posite cantilevered plates are: 

1. Laminate configurations can be selected for 
the cantilever that will result in the trailing 
edge tip deflecting either more or less than the 
leading edge tip where the load is applied. 

2. Failure to account for coupling due to non- 
symmetry can result in displacement differences 
that are of a magnitude of the order of the dis- 
placement. 

3. For the cases investigated, the vibration 
natural frequencies can be overestimated by as 
much as 90 percent if coupling due to nonsymmetry 
is not included, 

4. The vibration frequencies of the canti- 
levers investigated obtained using the 3 DOF per 
node idealization with reduced bending stiffness 
are within 3 to 7 percent of those obtained using 
the 5 DOF per node idealization. The displacements 
obtained by the same procedure are within 7 per- 
cent . 


5. Since the deformed shape of an unsymraetric 
cantilever when subjected to both loads and vibra- 
tions exhibits considerable bend-twist coupling, 
unsymmetric laminate configurations may be selected 
to yield predetermined deformation under load to 
provide built-in structural damping to minimize or 
avoid flutter. The deformed shape can be used to 
configure the laminate fabrication mold, 

6. The reduced stiffness approach using a 

5 DOF model as a baseline gives a good approxima- 
tion for determining displacements and vibration 
frequencies of imsyrametric anisotropic composite 
cantilevers . 

7. Finite element codes which do not account 
for bend-stretch coupling can be used .to predict 
displacements and vibration frequencies of unsym- 
metric laminates via the reduced bending stiffness 
approach ♦ 


Appendix - Derivation of Element Stiffness Matrix 

(Refer to Figs. 1 and 3) 

The displacement variables per node are u, v, 
w, ot, 3. Here u is taken along x, v along y, 

w along z; a is the rotation in the x-z plane 
and 3 in the y-z plane. The interpolation for- 
mula* 

6 

i=l 

defines a variable 0 within the element in terms 
of nodal variables 0^. The vector N is given by 






Li(2Li - 1) 
L2(2L2 - 1) 
L3(2L3 - 1) 
4 L 3 L 2 
4 L 2 L 3 


where the numbers refer to the similarly nuuhered 
nodes in Figiire 3 and where the area coordinates 
for parabolic variation are 



* 2^3 ■ ^ 3^2 ^2 ■ ^3 

- * 1^3 ^3 - 

^ 3 X 2 - * 2^1 - ^2 


24 = dET 


1 =^1 Xi 

1 X 2 72 
1 yn 


within the element the coordinates can be written 

6 


Svi 

i=l 


i=l 

6 




i=l 


Using the same Interpolation functions, the dis- 
placement variables are 


“ =2^ Vi 

i=l 

I 

6 

w = 

i=l 

6 

e =23 Vi 


o 

^=23 Vi 

i=l 

6 

“=23''i“i 


i-l 


i=l 


In terms of the area coordinates the strain- 
displacement relationships In familiar form are 
(subscripts denote direction. Fig. 3): 

e = e - zk 

XX XX XX 


yy 


e - zk 

yy yy 
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[K] 

30x30 


xy 


e - zk 
xy xy 


+ z^K 


xy 


where 


e = e 
yz yz 


3u 

3x 




yy 


5v 

3y 


w 

“2 


9u , 3v 

^ T — 

xy 3y Sx 


3x 


1 5u 
3x 


k « ig- - 

yy 3y 3y 


xy 


3y 3x 


1 3u ^ ^ 

Rl 3y R 2 3x 


xy “ ” R, 3y R^ 3x 


, 9w , u 
®xz = “ ^ 


yz 


« + in. + j!. 

® ^ ay ^ Rj 


The strain-displacement equations can be written 






^yy 

xy 


• [B] 

5x30 


5x1 


^6 

30x1 


The 8tresS“Strain equations can be written 

ro. 


XX 


< 0. 


"yy 

xy 


[E] 

5x5 


5x1 


^xx'* 


£yy 

^xy 

_^yz , 
5x1 


The element stiffness matrix can be written 


-I 


B^EB dV 


For a constant thickness thin plate the element 
stiffness matrix is 


[R] = 
30x30 


/ 


T * 

B D B dA 


where 



which are known as the plate constitutive equations. 
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Table 1 Summary of stiffnesses for a 12 -ply laminate with 

lT75-S/epoxy* FVR = 0,60] 



ious ply orientations 


(0,+45)3 

(-45,90)3 

(0,30)20,45 
-45 ,90 (-60, 90) 2 

(0,90)3 

«6 

(80, -40) 6 


VII 

viir 

IX 

X 

1783 

1872 

1783 

889 

556 

467 

556 

602 

0 

154 

553 

-618 

1783 

1872 

1783 

2585 

0 

-154 

553 

-94 

614 

525 

614 

660 

37 

37 

37 

37 

37 

37 

37 

37 

-32 

-43 

-18 

3.2 

0 

0 

13 

1.9 

-11 

-10 

14 

-2.7 

32 

43 

-8.7 

-7.0 

-11 

-10 

14 

-3.4 

0 

0 

13 

1.9 

0 

0 

0 


0 

0 

0 

0 

1597 

1679 

1716 

741 

352 

270 

463 

502 


142 

461 

-515 

1597 

1679 

1255 

2154 


-142 

461 

-78 


318 

511 

550 

31 

31 

31 

31 

31 

31 

31 

31 












Table 2 Summary of reduced bending stiffnesses 
for a 12-ply laminate with various ply 
orientations [T75-S/epoxy, FVR *» 0.60] 

(For ply orientation, see Table 1) 


Property 

Units 

Case 

IV 

VI 

VII 

VIII 

IX 

X 

I>11 

in. -lb 

374 

606 

733 

564 

682 

728 

^12 

in . -lb 

238 

6 

-66 

-52 

206 

492 

D13 

in. -lb 

0 

0 

-159 

-43 

202 

-504 

1^22 

In. -lb 

374 

606 

733 

564 

592 

2114 

I>23 

in. -lb 

0 

0 

159 

43 

189 

-77 

I>33 

in. -lb 

300 1 

68 

287 

231 

218 

540 

D44 

in. -lb 

31 

31 

31 i 

31 

31 

31 

^>55 

in . -lb 

31 

31 

31 

31 

31 

31 


Table 3 Tip deflection of a cantilevered unsymmetrlc fiber composite laminated plate (in.) 

(For ply orientation, see Table 1) 


(a) Without coupling (3 DOF, 16 element model) 


Location 

Case 




Aluminum 

I 

II 

■sa 

IV 

V 

VI 

VII 

VIII 

IX 

X 

LE 

0.0510 

0.0291 

0.0878 


0.0838 

0.0355 

0.0405 

0.0325 

0.0366 

0.0452 

0.108 

1/4 

,0483 

.0184 

.0869 

.0986 

.0837 

.0262 

.0306 

.2093 

.0319 

.0385 

,U6 

1/2 

.0453 

.0099 

.0851 


.0826 

.0172 

.0212 

.0264 

.0275 

.0326 

.124 

3/4 

,0432 

.0038 

.0825 

,0908 

.0807 

.0088 

.0123 

.0240 

.0234 

.0275 

.132 

TE 

,0408 

-,0014 : 

.0788 

.0856 

.0777 

.0007 

.0036 

.0217 

.0194 

.0229 

.140 



Finite element model 


(b) With coupling (5 DOF, 16 element model) 


Location 

Case 


Aluminum 

I 


III 

IV 

V 

VI 

VII 

VIII 

IX 

X 

LE 

0.0510 


mssi 




0.0887 

0.0562 

0,0771 

0.117 

B8H! S 

1/4 

.0483 

.0184 


.0986 

.146 

.0262 

.0762 

.0533 

.0724 

.101 

■Si 

1/2 

.0453 

,0099 

.0851 

.0950 

.143 


,0644 


.0685 

.087 


3/4 

.0432 

.0038 

.0825 


.139 

,0088 

.0536 

.0497 

.0656 

.075 


TE 

,0408 

-,0014 

.0788 

.0856 

.134 

.0007 

.0434 

.0487 

.0631 

.063 



(c) With reduced stiffness (3 DOF, 16 element model) 


Location 

Case 

IV 

VI 

VII 

VIII 

IX 

X 

LE 

0.146 

0,091 



0.109 


1/4 

.143 

.078 

,057 

.074 

,095 

.116 

1/2 

.139 

,066 

.055 

.070 

,082 

.124 

3/4 

.136 

.055 

.053 

.067 

,071 

.132 

TE 

.131 

.045 

.052 

.065 

.061 

.140 
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Table 4 Natural frequencies of a cantilevered unsynnDetrlc fiber composite 
laminated plate (cycles /sec) (For ply 
orientation see Table 1) 


(a) Without coupling (3 DOF, 16 -element model) 


Frequency 




Frequency ma^ltude for 

case 




order 

I 

II 

III 

IV 

V 

VI • 

VII 

VIII 

IX 

X 

1 

2 090 

824 

787 

831 

1 697 

1 535 

1 385 

1 407 

1 327 

606 

2 

2 859 

4 747 

4 524 

4 788 


2 562 

4 308 

3 947 

4 026 

3 207 

3 

8 521 

5 522 

5 076 

5 604 

8 810 

8 185 

7 567 

7 682 

7 346 

5 036 

4 

10 140 

12 304 

11 722 

12 411 

10 145 

9 724 

13 345 

12 476 

12 290 

8 804 

5 

U 022 

15 862 

14 422 

16 148 

20 355 

19 265 

18 091 

18 390 

17 774 

14 436 

6 


20 715 

20 200 

20 790 

20 960 

20 030 


22 370 


16 265 


(b) Without coupling (3 DOF, 4 -element model) 


Frequency 





case 




order 

I 

II 


IV 

V 


VII 

VIII 


X 



916 


936 

1 699 

1 537 

1 398 

1 411 

1 310 



2 786 

4 499 

4 025 

4 716 

2 565 

2 470 

4 160 

3 832 

4 069 



6 327 

5 329 

5 096 

5 300 

8 865 

8 259 

7 803 

7 636 

7 136 

5 125 


10 109 


10 186 

11 364 

11 879 

11 320' 

13 202 

13 248 

13 034 

9 780 


13 053 

15 077 

15 203 

14 754 

16 575 

16 203 


15 605 

15 626 

11 019 


14 738 

18 740 

17 084 


19 154 

20 398 

20 349 

20 371 

18 265 

21 141 


(c) With coupling (5 DOF, 4-element model) 


Frequency 

order 



VI 



IX 

X 

■M 

647 

892 

974 

867 


RQ 


3 152 

2 078 

3 445 

3 172 

2 867 



3 803 

5 322 

6 326 

5 551 

4 671 

sn 


7 207 

8 007 

10 438 

9 716 

9 601 

6 344 



8 359 

11 285 

11 296 

9 741 

9 727 

mBm 

10 811 

12 025 

14 374 

13 311 

12 124 

10 994 


Table 5 Comparison of the vibration frequencies of a cantilevered 
unsymmetric fiber composite laminated plate (case VII) with 
and without coupling and with reduced stiffness 


Frequency 


Frequency, cycles /sec 


order 

5 DOF /node 

3 DOF /node 

Bending only 


with coupling 

without coupling 

with coupling 


16 element 

4 element 

16 element 

4 element 

(approx) 

reduced stiffness 
4 element 

1 — 

960 

974 

1 385 

1 398 

931 


3 608 

3 445 

4 308 

4 160 

3 342 


5 611 

6 326 

7 567 

7 803 

6 128 


10 685 

10 438 

13 345 

13 202 

9 893 


11 335 

11 285 

18 091 

16 008 

13 811 


14 444 

14 374 

23 140 

20 349 

15 285 
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Table 6 Comparison of the vibration frequencies 
of a cantilevered unsymmetrlc fiber composite 
laminated plate with and without simulated 
centrifugal stiffening 

[Case VII, Table 1, 5 DOF, 16-element model] 


Frequency 

order 

Frequency , 

cycles /sec 

Percent 

difference 

Without 

stiffening 

With 

s tif f ening 

1 

974 

1 160 

20 

2 

3 445 

3 570 

4 

3 

6 326 

6 793 

7 

4 

10 438 

10 540 

1 

5 

11 285 

13 950 

24 

6 

14 374 

13 440 

7 


Table 7 Comparison of the tip displacements 
of a cantilevered unsymmetrlc fiber 
coniposite laminated plate 


[Case VII, Table I, 16-element model] 


Location 

Displacement, 

in. 

Without 
coupling 
(3 DOF) 

With 

coupling 
(5 DOF) 

Reduced 
stiffness 
(3 DOF) 

LE 

0,0325 

0,0562 

0.060 

1/4 

.0293 

.0533 

.057 

1/2 

,0264 

.0511 

.055 

3/4 

,0240 

,0497 

.053 

TE 

.0217 

.0487 

.052 
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Figure 1. - Schematic of the deformation (bending and twisting) of a 
cantilevered unsymmetric fiber composite laminated plate. 
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Figure 2. - Ply orientation schematic. 
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(a) DOUBLY CURVED VARIABLE (b) FLAT THIN PLATE (USED IN 
THICKNESS THIN SHELL PRESENT INVESTIGATION). 

Figure 3. - Schematic of the geometry of the 6-node triangular iso- 
parametric finite element. 
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a " 2 IN. 

(a) 4-ELEMENT MODEL (15-NODES). 

TRAILING EDGE 
MIDCHORD 
LEADING EDGE 

(b) 16-ElEMENT MODEL (45-NODES). 

Figure 4. - Finite element representation of cantilever. 





MODE ALUMINUM 
1 

2 

3 

4 

5 

6 



CASE I 


CASE VII 


CASE IX 

• 

1 

+ 

1 

+ 






1 - _ 


- 

1 



1 

+ 


+ 






$ + 

1 

1 

1 

- 1 ♦ 
f 

$ 

$ 

1 

- I* 







1 

) 

+ 

1 

+ ; ^ 

§ - 1 + 
§ / 

1 

§ 

" ' - T* 

1 

'.V. 

- / ♦ 





-i-. 

-O 1 

1 

1 



- / / ^ 

y 1 

♦ / - ; + 

1. - 

1 

- '\ * 

1 







1— =— 

1 

* ! - '<* 
t j 

1 

,N 

1 

/ 

d ljL 


Figure 5. - First six mode shapes of a cantilevered unsymmetric fiber 
composite laminated plate. (For case identification, see Table I. ) 
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